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Introduction:  Prostate  cancer  (PCa)  is  the  most  common  cancer  in  men  in  several 
countries,  with  the  American  Cancer  Society  (ACS)  estimating  241,740  new  cases  of  PCa 
to  be  diagnosed  and  deaths  of  28,170  men  of  PCa  (1).  Due  to  its  prevalence  in  the  male 
population  as  well  as  its  unpredictable  clinical  course,  early  detection  and  diagnosis  have 
become  a  priority  for  many  health  care  professionals.  Another  method  for  staging 
prostate  cancer  is  through  imaging  techniques  including  ultrasound,  computed 
tomography  (CT),  and  magnetic  resonance  imaging  (MRI)  with  or  without  the  help  of 
dynamic  contrast  enhancement  modeling  (DCE-MRI),  diffusion  weighted  imaging 
(DWI),  and  magnetic  resonance  spectroscopy  (MRS)  (2-5).  MRS  is  a  powerful  tool  for 
exploring  the  cellular  chemistry  of  human  tissues  (3,5,6-11).  There  is  a  growing  body  of 
evidence  that  1H  MRS  may  contribute  to  the  clinical  evaluation  of  prostate  cancer  and 
also  for  evaluating  the  metabolic  alterations  due  to  therapy.  There  have  been  no  reports 
on  combining  two  spectral  dimensions  with  two-dimensional  (2D)  or  three  dimensional 
(3D)  spatial  encoding  applicable  to  prostate  cancer.  Acceleration  of  magnetic  resonance 
spectroscopic  imaging  (MRSI)  has  been  demonstrated  using  echo-planar  imaging 
techniques  (12-13).  Recently,  Schulte  et  al.  have  successfully  developed  an  algorithm 
called  prior-knowledge  fitting  (ProFit)  to  quantify  metabolite  concentrations  using  the 
JPRESS  spectra  recorded  using  a  Philips  3T  MRI  scanner  (14).  It  was  demonstrated  that 
metabolite  quantitation  of  JPRESS  spectra  with  ProFit  was  accurate,  robust  and  yielding 
generally  consistent  results,  both  in  vivo  and  in  vitro.  Their  results  suggest  that  the 
number  of  quantifiable  prostate  metabolites  can  be  increased  from  3-4  with  ID 
PRESS/LC-Model  to  more  than  10  with  JPRESS/ProFit  (15-16). 

Body: 

i)  Proposed  Task  1  (Months  1-6):  To  implement  a  multi-voxel  based  extension  of  the 
JPRESS  sequence,  in  which  two  spectral  encodings  will  be  combined  with  two  spatial 
encodings  using  the  new  Siemens  VB17a  platform.  This  four-dimensioned  (4D)  data 
acquisition  scheme  will  be  accomplished  utilizing  the  EPI  approach  that  is  commonly 
used  for  spatied  encoding  in  MRI. 

Completed  and  Reported  in  the  1st  year  Annual  Report. 

ii)  Proposed  Task  2:  To  eveduate  the  EPI-bcised  JPRESS  using  a  prostate  phantom 
containing  severed  meteibolites  which  have  been  reported  in  prostate  tissues,  and  to 
optimize  the  EP-JJRESI  sequence  anel  other  acquisition  parameters  using  the  phantom 

(Months  6-12). 

Completed  and  Reported  in  the  1st  year  Annual  Report. 

iv)  Proposed  Task  3:  To  record  the  4D  EP-JRESI  spectra  in  the  peripheral,  centred  eind 
treinsition  zones  of  healthy  prostates.  (Months  6-18). 

Five  more  healthy  males  have  been  investigated. 

Proposed  Task  4:  To  elevelop,  eveduate  eind  optimize  the  prior -knowleelge  basis  set 
spectra  using  the  GAMMA- simulation  anel  prostate  phantom  solutions  as  prior 
knowleelge  for  the  multi-voxel  based  JPRESS  spectra  recoreled  using  the  3T  MRI  scanner 

(Months  6-12). 
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Accomplished  during  September  29,  2011-October  28  2012: 


Proposed  Task  5:  To  record  the  multi- voxel-based  2D  spectra  in  the  peripheral,  central 
and  transition  zones  of  patients  with  BPH  and  malignant  prostate  cancer.  The  prostate 
metabolite  concentrations  calculated  using  the  ProFit  algorithm  prepared  for  the  multi¬ 
voxel  data  will  be  compared  with  LC-Model  processed  ID  spectral  based  MRSI  data 

(Months  18-48). 


Accomplished  during  October  29,  2014-October  28  2015:  Ten  malignant  prostate 
cancer  patients  and  three  healthy  males  were  investigated  during  the  current  year.  The  4D 
EP-JRESI  data  were  recorded  using  the  following  parameters:  TR/TE=1.5s/30ms,  2 
averages,  512  ti,  oversampled  32kx,  64  increments  along  the  indirect  spectral  (ti)  and  16 
spatial  ky  dimensions;  the  endorectal  "receive"  coil.  The  oversampled  32kx  were 
processed  to  normal  16kx  points.  First,  the  4D  NUS  EP-JRESI  data  was  acquired  in  a  71 
year  old  PCa  patient;  In  our  recent  publication,  we  had  included  twenty  two  PCa  patients 
with  a  mean  age  of  63.8  years  (range:  46-79  years)  who  subsequently  underwent  radical 
prostatectomy.  Patients’  Gleason  scores  varied  between  6  and  9.  Prostate-specific  antigen 


(A) 


(B) 


Figure  1.  Spatial  maps  of  (Ch  +Cr)  for  TV  (A)  and  MaxEnt  (B)  non-linear  reconstruction 
methods  of  the  4D  EP-JRESI  data  recorded  in  a  74  year  old  PCa  patient.  The  EP-JRESI 
was  overlaid  on  top  of  the  T2  weighted  MRI 

We  were  able  to  reproduce  the  above  mentioned  findings  in  more  patients  with 
known  GS,  and  also,  in  remaining  four  patients  whose  histology  findings  were  not 
known.  2D  peak  (cross  and  diagonal)  volumes  were  defined  by  the  operator  in  the 
frequency  domain  2D  JPRESS  data  for  this  pilot  analysis.  Two  different  non-linear 
reconstruction  methods  such  as  maximum  entropy  (MaxEnt)  and  total  variation  (TV).  As 
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evident  in  Fig. 2,  selected  ratios  showed  significant  changes  (p<0.05);  however,  this 
outcome  from  a  small  group  of  PCa  patients 

Fig-2.  Metabolite  ratios  of  Cit,  Spm,  ml  and  Glx  in  cancer  and  non-cancer  locations 
processed  by  non-linear  reconstruction  methods  using  TV  and  MaxEnt. 
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needs  to  be  treated  with  caution  showing  a  necessity  for  further  research  using  larger 
patient  cohorts. 


Key  Research  Accomplishments 


•  Continued  evaluation  of  the  4D  EP-JRESI  sequence  in  10  malignant  prostate 
cancer  patients  and  three  healthy  male  subjects. 

•  We  have  made  further  progress  on  compressed  sensing  reconstruction  of  the 
non-uniformly  undersampled  4D  EP-JRESI  sequence  using  different 
reconstruction  methods  such  as  maximum  entropy,  total  variation  (TV),  etc  .  This 
will  facilitate  shortening  the  endorectal  spectral  acquisition  and  reducing  the 
patient  inconvenience  during  the  scan.  A  peer-reviewed  manuscript  was  published 
recently  in  the  journal  of  ‘NMR  in  Biomedicine’. 

•  Using  the  preliminary  results  obtained  using  this  IDEA  grant,  we  had  submitted  a 
NIH  R01  grant  application  submission  entitled  “Fast  J-resolved  Prostate  MR 
Spectroscopic  Imaging  and  Non-linear  Reconstruction”  was  resubmitted  on 
February  5,  2015  and  we  plan  to  resubmit  this  grant  to  NIH. 

Reportable  Outcomes: 

A.  Peer-reviewed  Publications: 

1.  Nagarajan  R,  Iqbal  Z,  Burns  B,  et  al.  Accelerated  Echo  Planar  J-Resolved 

Spectroscopic  Imaging  in  Prostate  Cancer:  Nonlinear  Reconstruction  Using  Total 

Variation  And  Maximum  Entropy.  NMR  in  Biomedicine  2015;  28:  1366-1373. 


B.  Presentations:  During  the  4th  year,  two  abstracts  were  presented:  (i)  first  at  the  23ld 
International  Society  of  Magnetic  Resonance  in  Medicine  (ISMRM),  Toronto,  Canada, 
May  30-June  5,  2015:  the  abstract  entitled  “Semi-Laser  5D  Echo-Planar  J-Resolved 
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Spectroscopic  Imaging:  Pilot  Validation  in  Prostate  Cancer.”  (ePoster  presentation), 
#2020;  the  2nd  abstract  entitled  “How  Accurately  Multiparametric  MRI  Detect  Prostate 
Cancer?”  (eposter),  #  LL-PHE1 18;. 

C.  Books:  A  book  chapter  including  some  prostate  data  was  published  in  the  following 
chapter:  Huda  A,  Nagarajan  R,  Furuyama  J  and  Thomas  MA.  In  Vivo  Two-Dimensional 
Magnetic  Resonance  Spectroscopy  "Comprehensive  Biomedical  Physics"  Elsevier 
Publishers,  2014,  Pp  359-77. 

Conclusions:  After  successfully  implementing  the  4D  EP-JRESI  scanning  protocol  on 
the  3T  MRI  scanner  and  testing  it  in  healthy  controls  and  patients  during  the  earlier  3 
years,  13  more  subjects  have  been  investigated  in  the  4th  year.  We  will  continue  to 
recruit  10  malignant  and  benign  prostate  cancer  patients,  and  5  healthy  male  subjects 
during  the  next  year  (2nd  no-cost  time  extension).  Two  new  manuscripts  and  few 
conference  abstracts  will  be  prepared  to  present  the  summary  of  findings. 
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Accelerated  echo  planar  ./-resolved 
spectroscopic  imaging  in  prostate  cancer:  a 
pilot  validation  of  non-linear  reconstruction 
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The  overlap  of  metabolites  is  a  major  limitation  in  one-dimensional  (ID)  spectral-based  single-voxel  MRS  and 
multivoxel-based  MRSI.  By  combining  echo  planar  spectroscopic  imaging  (EPSI)  with  a  two-dimensional  (2D) 
./-resolved  spectroscopic  (JPRESS)  sequence,  2D  spectra  can  be  recorded  in  multiple  locations  in  a  single  slice  of  pros¬ 
tate  using  four-dimensional  (4D)  echo  planar  ./-resolved  spectroscopic  imaging  (EP-JRESI).  The  goal  of  the  present 
work  was  to  validate  two  different  non-linear  reconstruction  methods  independently  using  compressed  sensing- 
based  4D  EP-JRESI  in  prostate  cancer  (PCa):  maximum  entropy  (MaxEnt)  and  total  variation  (TV).  Twenty-two 
patients  with  PCa  with  a  mean  age  of  63.8  years  (range,  46-79  years)  were  investigated  in  this  study.  A  4D 
non-uniformly  undersampled  (NUS)  EP-JRESI  sequence  was  implemented  on  a  Siemens  3-T  MRI  scanner.  The  NUS 
data  were  reconstructed  using  two  non-linear  reconstruction  methods,  namely  MaxEnt  and  TV.  Using  both  TV  and 
MaxEnt  reconstruction  methods,  the  following  observations  were  made  in  cancerous  compared  with  non- 
cancerous  locations:  (i)  higher  mean  (choline  +  creatine)/citrate  metabolite  ratios;  (ii)  increased  levels  of  (choline  + 
creatine)/spermine  and  (choline  +  creatine)/myo-inositol;  and  (iii)  decreased  levels  of  (choline  +  creatine)/(glutamine 
+  glutamate).  We  have  shown  that  it  is  possible  to  accelerate  the  4D  EP-JRESI  sequence  by  four  times  and  that  the 
data  can  be  reliably  reconstructed  using  the  TV  and  MaxEnt  methods.  The  total  acquisition  duration  was  less  than 
13  min  and  we  were  able  to  detect  and  quantify  several  metabolites.  Copyright  ©  2015  John  Wiley  &  Sons,  Ltd. 
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INTRODUCTION 

Prostate  cancer  (PCa)  is  the  most  commonly  diagnosed  non- 
cutaneous  malignancy  in  the  USA  and  is  the  second  leading 
cause  of  cancer-related  death  in  men  (1).  One  in  six  men  will 
be  diagnosed  with  PCa  during  their  lifetime,  but  only  one  in  36 
will  die  of  this  disease.  Currently,  the  annual  prostate-specific  an¬ 
tigen  (PSA)  test  and  digital  rectal  examination  (DRE)  are  routinely 
performed  (2)  for  screening.  The  PSA  screening  test  measures 
the  serum  level  of  PSA  in  blood  samples.  However,  it  is  a  contro¬ 
versial  test  because  65-75%  of  PSA  screening  gives  false-positive 
results  leading  to  overdiagnosis  (3).  The  use  of  systematic 
transrectal  biopsy  can  miss  significant  cancer  lesions  because 
of  random  sampling  error  (4)  and  the  observation  that  one-third 
of  significant  tumors  lie  in  the  anterior  part  of  the  gland,  based 
on  studies  of  radical  prostatectomy  specimens  (5).  Hence,  there 
is  an  immediate  need  for  early,  yet  accurate,  detection  of  PCa 
to  improve  disease  outcomes. 

nH  MRS  enables  the  detection  of  a  range  of  biochemicals  in 
the  prostate  by  making  use  of  the  proton  signals  in  these 
molecules.  The  detection  of  biochemicals  in  vivo  is  limited  to 
concentrations  of  more  than  0.5-1  rriM.  Signals  of  citrate  (Cit), 
creatine  (Cr),  choline  (Ch)  and  spermine  (Spm)  can  be  detected 


throughout  the  prostate,  with  increased  levels  of  Ch  and  de¬ 
creased  levels  of  Cit  being  indicative  of  cancer  (6-8). 
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Current  limitations  of  single-voxel-based  MRS  and  MRSI  in  the 
prostate  are  caused  by  the  overlap  of  metabolite  resonances, 
allowing  the  quantification  of  only  a  few  metabolites  (Cit,  Ch, 
Cr  and  Spm)  and  the  use  of  long  TEs.  The  conventional  MRSI 
technique  can  be  accelerated  by  echo  planar  spectroscopic 
imaging  (EPSI)  (9-13).  EPSI  speeds  up  MRSI  using  an  echo  planar 
readout  of  one  spectral  and  one  spatial  dimension,  thereby 
achieving  an  acceleration  factor  equal  to  the  number  of  points 
along  one  of  the  spatial  dimensions.  For  example,  a  two- 
dimensional  (2D)  spatial  matrix  array  (16  x  16)  would  be 
acquired  16  times  faster  with  EPSI  than  with  conventional  MRSI. 
Flowever,  the  acceleration  may  be  at  the  cost  of  the  signal-to- 
noise  ratio  (SNR)  (11)  and  the  spectra  could  be  affected  by 
Nyquist  ghost  artifacts  (14). 

A  single-voxel-based  2D  J-resolved  spectroscopic  (JPRESS) 
sequence  has  been  evaluated  in  PCa,  and  has  shown  improved 
spectral  dispersion  because  of  the  added  spectral  dimension 
(15,16).  New  computational  methods  have  made  compressed 
sensing  feasible  to  accelerate  MRI  by  exploiting  the  sparsity  of 
the  images  in  a  known  transform  domain  to  reconstruct  non- 
uniformly  undersampled  (NUS)  k-space  data  (17).  For  further 
acceleration,  the  application  of  compressed  sensing  for  MRSI  is 
apt,  exploiting  sparsity  in  multiple  dimensions  of  frequency 
and  space  in  transform  domains  of  wavelets  and  total  variation 
(TV)  (18-20).  By  combining  EPSI  with  JPRESS,  2D  spectra  can  be 
recorded  in  multiple  locations  in  the  prostate  using  four¬ 
dimensional  (4D)  echo  planar  J-resolved  spectroscopic  imaging 
(EP-JRESI),  which  combines  two  spectral  with  two  spatial  dimen¬ 
sions.  A  pilot  feasibility  has  been  demonstrated  recently  to  map 
metabolites  in  the  healthy  human  prostate  and  brain  (21,22). 

Maximum  entropy  (MaxEnt)  reconstruction  finds  the  spectrum 
that  maximizes  entropy  whilst  maintaining  consistency  with  the 
measured  data.  MaxEnt  reconstruction  is  an  alternative  non-linear 
reconstruction  technique  to  compressed  sensing.  MaxEnt  has 
been  successfully  used  to  reconstruct  undersampled  images  in 
astronomy  and  multidimensional  spectra  in  NMR  (23-25),  but  has 
not  been  applied  to  the  spatial-spectral  domain  (ky-fq)  of  4D 
EP-JRESI  of  PCa.  MaxEnt  and  TV  algorithms  have  been  used  to  re¬ 
construct  the  NUS  indirect  spectral  and  spatial  dimensions  (21 ,26). 

The  TV  algorithm  was  first  proposed  by  Rudin  ef  at.  (27)  for  im¬ 
age  denoising  and,  since  then,  has  been  successfully  used  for  im¬ 
age  restoration.  In  the  TV  algorithm,  an  objective  function  using 
the  TV  norm  is  minimized  subject  to  a  data  fidelity  term  posed 
by  the  acquired  projection  data.  Minimization  of  the  image  gradi¬ 
ent  essentially  suppresses  those  high  spatial  frequency  parts, 
such  as  streaking  artifacts  and  noise,  in  the  reconstructed  images. 

The  goal  of  the  present  work  was  to  validate  the  MaxEnt  and 
TV  non-linear  reconstruction  algorithms  separately  in  patients 
with  PCa  using  compressed  sensing-based  4D  EP-JRESI  data. 

MATERIALS  AND  METHODS 

Patients 

Between  March  2012  and  May  2013,  twenty-two  patients  with 
PCa  with  a  mean  age  of  63.8  years  (range,  46-79  years),  who  sub¬ 
sequently  underwent  radical  prostatectomy,  were  selected  for 
the  study.  The  patients'  Gleason  scores  varied  between  6  and 
9.  Their  PSA  levels  varied  from  0.7  to  22.8  ng/mL  (mean,  6.23 
ng/mL).  These  patients  were  scanned  using  a  3-T  Siemens 
(Siemens  Medical  Solutions,  Erlangen,  Germany)  MRI  scanner 
with  an  endorectal  'receive'  coil.  The  protocol  combining  MRI 


and  MRS  was  performed  at  least  8  weeks  after  transrectal 
ultrasound-guided  sextant  biopsy.  The  entire  protocol  was  ap¬ 
proved  by  the  Institutional  Review  Board,  and  informed  consent 
was  obtained  from  each  patient.  PCa  was  histopathologically 
confirmed  after  radical  prostatectomy.  The  voxels  covering  the 
tumorous  lesions  from  the  peripheral  zone  (PZ)  were  selected 
and  indicated  as  tumor  voxels,  which  was  confirmed  by  the 
pathology  report.  After  reconstruction,  the  EP-JRESI  data  were 
overlaid  onto  MRI  images. 

MRI  and  MRSI 

A  body  matrix  phased-array  coil  assembly,  combined  with  an 
endorectal  coil,  was  used  in  the  'receive'  mode,  whereas  a  quad¬ 
rature  body  'transmit'  coil  was  used.  All  patients  were  imaged  in 
the  supine  (feet-first)  position.  Axial  images  were  oriented  to  be 
perpendicular  to  the  long  axis  of  the  prostate,  which  was  guided 
by  the  sagittal  images.  Axial,  coronal  and  sagittal  T2-weighted 
(T2W)  turbo  spin  echo  images  were  recorded  using  the  following 
parameters:  TR/TE  =  3850-4200  ms/96-101  ms;  slice  thickness, 
3  mm;  field  of  view,  20  x  20  cm2;  echo  train  length,  13;  data 
matrix,  320  x  256. 

A  compressed  sensing-based  4D  EP-JRESI  sequence  was  vali¬ 
dated  on  the  3-T  MRI  scanner  and  the  volume  of  interest  (VOI) 
was  localized  using  three  slice-selective  radiofrequency  (RF) 
pulses  (90°-180°-180°)  (Fig.  1).  The  total  time  for  the  acquisition 
of  a  fully  sampled  4D  EP-JRESI  scan  (TR  =  1.5  s,  16ky  x  16kx, 
64-1 00fq,  512 12)  can  be  more  than  25  min.  The  parameters  for 
EP-JRESI  were  as  follows:  TR/TE/Avg  =  1500  ms/30  ms/2;  16 
phase-encoding  steps;  512  complex  points  with  an  F2  bandwidth 
of  1 1 90  Hz  along  the  detected  dimension.  For  the  second  dimen¬ 
sion  (Ft ),  64  increments  with  bandwidths  of  1000  Hz  were  used. 
The  in-plane  spatial  resolution  and  slice  thickness  were  1  x  1  cm2 
and  1  cm,  respectively.  As  the  EPSI  readout  simultaneously  ac¬ 
quires  one  spatially  encoded  dimension  (kx)  and  one  temporal 
dimension  (f2),  we  propose  the  use  of  NUS  in  the  remaining 
ky-f,  plane,  followed  by  compressed  sensing  reconstruction 
(MaxEnt  and  TV).  A  4  x  NUS  was  imposed  along  the  plane 
containing  the  incremented  spectral  and  spatial  dimensions 
and  ky).  Despite  the  mixing  of  spatial  and  spectral  dimensions 
in  the  reconstruction,  the  sparsity  requirement  for  reconstruction 
is  shown  to  be  satisfied,  as  required  by  compressed  sensing.  As 
the  ky-ti  values  are  incremented,  NUS  can  be  applied  along 
the  ky-ti  plane. 

The  individual  voxel  volume  in  human  prostate  was  1  mL.  Two 
sets  of  data  were  collected:  one  water-suppressed  (WS)  scan  with 
a  total  scan  time  of  12  min  and  a  second  non-water-suppressed 
(NWS)  scan  using  one  average  and  one  f,  increment  (30  s).  The 
NWS  scan  was  used  for  phase  corrections  (eddy  current  correc¬ 
tions).  The  full  width  at  half-maximum  of  the  water  peak  in  the  can¬ 
cerous  and  non-cancerous  locations  was  between  20  and  25  Hz. 

Data  analysis 

The  NUS  data  were  reconstructed  by  MaxEnt  and  TV  separately. 
A  modified  Split-Bregman  algorithm  (28)  solves  the  uncon¬ 
strained  TV  optimization  problem  as: 

min  ||  Vm  ||  t  +  2||Fum  —  y||2 

m 

where  V  is  the  gradient  operator,  m  is  the  reconstructed  data, 
||x|h  is  the  /q-norm,  X  is  a  regularization  parameter,  Fu  is  the 
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Figure  1.  Schematic  diagram  of  the  four-dimensional  (4D)  echo  planar  ./-resolved  spectroscopic  imaging  (EP-JRESI)  pulse  sequence  with  the  volume  of 
interest  (VOI)  localized  using  point-resolved  spectroscopy  with  three  radiofrequency  (RF)  pulses  (90°,  180°,  180°).  The  indirect  dimension  (fn)  was  intro¬ 
duced  before  the  last  180°  pulse.  ADC,  analog  to  digital  converter;  Gs,  GP  and  GF,  gradients  for  slice  selection,  phase  encoding  and  frequency  encoding, 
respectively;  f2,  detected  dimension;  n,  number  of  echo  planar  imaging  bipolar  readouts. 


undersampled  Fourier  transform  and  y  is  the  undersampled 
data.  The  above  equation  removes  the  incoherent  artifacts 
caused  by  NUS  by  minimizing  TV,  whilst  maintaining  consis¬ 
tency  with  the  sampled  measurements.  The  TV  regularization 
parameters  were  the  same  as  reported  by  Burns  et  at.  (29). 
We  have  used  the  Split-Bregman  reconstruction  method  pri¬ 
marily  for  its  robustness  against  the  regularization  parameters 
chosen.  Because  of  the  use  of  Bregman  parameters  in  the 
reconstruction  algorithm,  which  are  calculated  using  the  differ¬ 
ence  between  the  reconstruction  and  the  sampled  data  at 
each  iteration,  the  influence  of  the  regularization  parameters 
is  greatly  lessened  compared  with  the  use  of  other  algorithms 
that  solve  the  TV  problems.  Although  the  choice  of  parameters 
can  influence  the  reconstruction,  the  algorithm  allows  for  a 
wider  range  of  possible  values  in  order  to  achieve  roughly 
the  same  results. 

MaxEnt  is  a  constrained  convex  optimization  algorithm  that 
uses  a  variant  of  the  conjugate  gradient  method  to  iteratively 
solve  the  inverse  problem  (23,25,28): 

maximize  Sq/2(f)  such  that  || F_1  Kf  —  D\\2<a  [2] 

where  f  is  the  estimated  fully  sampled  spectrum  at  each  itera¬ 
tion,  F1  is  the  inverse  Fourier  transform,  K  is  the  NUS  matrix,  D 
is  the  measured  time  domain  data,  a  is  the  noise  standard 
deviation  and  SU2(f)  is  the  spin  -  Vi  entropy  of  the  estimated 
spectrum  (24).  All  compressed  sensing  4D  EP-JRESI  data  were 
processed  using  TV  and  MaxEnt  reconstruction  with  custom 
MATLAB  software.  The  reconstruction  time  for  each  method  took 
about  25  min  using  an  8  GB  RAM,  Intel  Core  i7-2600  CPU  @  3.40 
GFIz.  For  the  2D  data  processing,  the  raw  matrix  was  apodized 
with  phase  shifted  and  squared  sine  bell  functions  along  fq  and 
f2,  and  zero  filled  to  128  x  1024  prior  to  fast  Fourier  transforma¬ 
tion  along  the  two  dimensions.  All  2D  spectra  were  presented  as 
contour  plots,  and  the  2D  spectral  matrices  were  not  skewed  by 
45°  about  Fq  =  0  FHz. 

Statistical  analysis 

Statistical  analyses  were  performed  with  SPSS  21  (SPSS  Inc., 
Chicago,  IL,  USA).  Using  logistic  regression  analysis,  areas  under 
the  curve  (AUCs)  of  the  receiver  operating  characteristic  (ROC) 
were  calculated  for  various  metabolites  to  discriminate  between 


MaxEnt  and  TV  reconstruction  methods.  In  addition,  the  paired 
f-test  was  used  to  determine  the  various  metabolite  ratios  in 
cancerous  and  non-cancerous  locations,  p  <  0.05  was  consid¬ 
ered  to  be  statistically  significant. 


RESULTS 

Using  this  pilot  validation,  2D  peaks  attributed  to  Cit,  Ch,  Cr,  Spm, 
myo-inositol  (ml)  and  glutamate  (Glu)  plus  glutamine  (Gin)  (Glu  + 
Gin  =  Glx)  were  quantified  in  cancerous  and  non-cancerous  loca¬ 
tions  using  the  peak  integration  MATLAB  code.  Figure  2  shows 
(Ch  +  Cr)/Cit,  (Ch  +  Cr)/Spm,  (Ch  +  Cr)/ml  and  (Ch  +  Cr)/Glx  of 
cancerous  and  non-cancerous  locations  processed  by  TV  and 
MaxEnt.  The  mean  metabolite  ratios  (±  standard  deviation,  SD) 
of  Cit,  Spm,  ml  and  Glx  of  the  non-cancerous  locations,  proc¬ 
essed  using  TV,  were  1.158  ±  0.830,  2.396  ±  1.95,  5.325  ±  2.42 
and  5.404  ±  2.74,  respectively.  In  the  cancerous  locations,  the  cor¬ 
responding  metabolite  ratios  were:  4.209  ±  2.132,  2.808  ±  1.77, 
5.640  ±  2.18  and  5.275  ±  2.80.  Similarly,  the  mean  metabolite 
ratios  (±SD)  of  Cit,  Spm,  ml  and  Glx  of  the  non-cancerous  loca¬ 
tions,  calculated  using  the  MaxEnt-reconstructed  data,  were 
1.079  ±  0.795,  2.096  ±  1.06,  4.967  ±  2.114  and  5.902  ±  3.40, 
respectively.  In  the  cancerous  locations,  the  corresponding  me¬ 
tabolite  ratios  were  3.620  ±  1.759,  2.727  ±  1.46,  6.008  ±  2.57 
and  5.275  ±  3.19. 

We  found  that  the  mean  Cit  metabolite  ratios  were  signifi¬ 
cantly  higher  in  cancerous  locations  relative  to  non-cancerous  lo¬ 
cations  in  both  the  TV  and  MaxEnt  reconstructions  (p  <  0.005). 
Increased  levels  of  Spm  (p  =  0.46)  and  ml  (p  =  0.65)  ratios,  and 
decreased  levels  of  Glx  (p  =  0.88)  ratios,  were  observed  in  cancer¬ 
ous  locations  relative  to  non-cancerous  locations  in  the  TV  recon¬ 
struction.  Similarly,  in  the  MaxEnt  reconstruction,  increased 
levels  of  Spm  (p  =  0.25)  and  ml  (p  =  0.15)  ratios,  and  decreased 
levels  of  Glx  (p  =  0.81)  ratios,  were  observed  in  cancerous 
locations  relative  to  non-cancerous  locations.  None  of  the  ratios 
could  discriminate  significantly  between  differing  grades 
(Gleason  scores)  of  PCa  because  of  overlap  of  the  ratio  values. 

Figure  3  shows  spatial  maps  of  (Ch  +  Cr)  for  the  TV  and  MaxEnt 
reconstructed  data  acquired  in  a  74-year-old  patient  with  PCa. 
The  MaxEnt  and  TV  reconstructions  of  the  cancerous  (Fig.  4B, 
D)  and  non-cancerous  (Fig.  4C,  E)  locations  extracted  from  the 
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Figure  2.  Metabolite  ratios  of  citrate  (Cit),  spermine  (Spm),  myo-inositol  (ml)  and  glutamine  +  glutamate  (Glx)  in  cancerous  and  non-cancerous  loca¬ 
tions  processed  by  non-linear  reconstruction  methods  using  total  variation  (TV)  and  maximum  entropy  (MaxEnt).  Ch,  choline;  Cr,  creatine. 


JPRESS  spectrum,  obtained  from  the  4D  EP-JRESI  data,  are  shown 
in  Fig.  4.  Figure  4C  illustrates  the  regions  of  interest  used  for  peak 
integration.  We  compared  and  correlated  the  TV  and  MaxEnt 
reconstruction  methods  for  Cit,  Spm,  ml  and  Glx  in  the  cancerous 
and  non-cancerous  locations.  The  correlation  of  the  (Ch  +  Cr)/Cit, 
(Ch  +  Cr)/Spm,  (Ch  +  Cr)/ml  and  (Ch  +  Cr)/Glx  ratios  for  the 
MaxEnt  and  TV  reconstructions  in  the  cancerous  locations  are 
shown  in  Fig.  5.  For  each  patient,  two  to  three  voxels  were 
selected  in  the  PZ  of  the  cancerous  and  non-cancerous  locations, 
and  the  average  values  for  each  location  were  reported.  A  posi¬ 
tive  correlation  was  found  for  the  following  metabolites  in  the 
cancerous  locations:  (Ch  +  Cr)/Cit  (R2  =  0.85),  (Ch  +  Cr)/Glx  ( R 2 
=  0.96),  (Ch  +  Cr)/Spm  (R2  =  0.86)  and  (Ch  +  Cr)/ml  (R2  =  0.95). 
The  concentration  of  Cit  is  higher  in  healthy  prostate.  Flence,  if 
the  Cit  peak  was  higher  than  the  Ch  peak,  the  voxel  was  consid¬ 
ered  to  be  non-cancerous  for  (Ch  +  Cr)/Cit  values  below  0.5  and 
malignant  for  (Ch  +  Cr)/Cit  values  above  0.5.  These  values  were 
selected  manually  on  each  subject. 

The  results  of  the  logistic  regression  analysis  and  consequent 
ROC  curve  analyses  are  given  in  Table  1,  including  the  sensitivity, 


specificity,  positive  predictive  value  (PPV),  negative  predictive 
value  (NPV),  AUC  and  accuracy  for  the  classification  of  the 
MaxEnt  and  TV  methods.  ROC  curve  analyses  for  differentiating 
the  metabolite  ratios  of  cancerous  and  non-cancerous  locations 
in  MaxEnt  suggest  that  the  Cit  ratio  gives  the  best  predictability, 
with  a  sensitivity  of  86.4%,  specificity  of  90.0%,  accuracy  of  88.6% 
and  AUC  =  94.0%.  In  addition,  the  overall  sensitivity,  specificity, 
accuracy  and  AUC  of  ml  and  Spm  were  slightly  better  in  MaxEnt 
compared  with  TV. 


DISCUSSION 

Using  the  NUS  data  with  non-linear  iterative  reconstruction,  we 
have  validated  the  TV  and  MaxEnt  reconstruction  methods 
independently  in  patients  with  PCa  using  EP-JRESI  in  a  clinically 
feasible  time.  In  addition  to  a  significantly  increased  Cit  ratio  in 
cancerous  locations,  increased  metabolite  ratios  of  Spm  and  ml, 
and  decreased  ratios  of  Glx,  were  found  in  cancerous  locations 
compared  with  non-cancerous  locations.  Although  TV  and 


Figure  3.  Spatial  maps  of  (choline  +  creatine)  [(Ch  +  Cr)]  for  total  variation  (TV)  (A)  and  maximum  entropy  (MaxEnt)  (B)  non-linear  reconstruction 
methods  of  the  four-dimensional  (4D)  echo  planar  J-resolved  spectroscopic  imaging  (EP-JRESI)  data  recorded  in  a  74-year-old  patient  with  prostate 
cancer  (PCa).  EP-JRESI  was  overlaid  on  top  of  the  T2-weighted  MRI. 
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Figure  4.  Echo  planar  ./-resolved  spectroscopic  imaging  (EP-JRESI)  voxel  localization  on  top  of  the  T2~weighted  MRI  (A);  two-dimensional  (2D) 
./-resolved  spectroscopy  (JPRESS)  spectra  extracted  from  the  maximum  entropy  (MaxEnt)  and  total  variation  (TV)  reconstructions  of  cancerous  (B,  D)  and 
non-cancerous  (C,  E)  locations.  The  regions  of  interest  used  for  peak  integration  are  shown  in  (C);  (MM,  macromolecules;  Cit,  citrate;  Glx,  glutamine  + 
glutamate;  ml/PCh,  myo-inositol/phosphocholine;  Tau,  taurine;  Spm,  spermine). 


MaxEnt  reconstruction  methods  showed  comparable  results  in 
cancerous  and  non-cancerous  locations,  the  sensitivity,  accuracy 
and  AUC  were  slightly  increased  in  the  MaxEnt  reconstruction. 

In  the  present  study,  we  report  the  ratios  of  (Ch  +  Cr)/Cit, 
(Ch  +  Cr)/Spm,  (Ch  +  Cr)/ml  and  (Ch  +  Cr)/Glx,  because  of  the 
proximity  of  the  total  Cr  peak  (3.0  ppm)  to  the  total  Ch  peak 
(3.2  ppm)  in  these  in  vivo  MR  spectra,  which  were  therefore 
not  always  separable.  In  this  study,  significantly  higher  ratios 
of  Cit  were  observed  in  the  PZ  of  cancerous  locations. 

It  is  likely  that  the  drop  in  Cit  levels  precedes  malignant  trans¬ 
formation  (30).  It  has  been  suggested  that,  as  a  result  of  a  meta¬ 
bolic  switch,  neoplastic  cells  oxidize  Cit,  whereas  normal 
prostatic  cells  show  a  low  Cit  oxidizing  capability  (31).  Decreased 
levels  of  zinc,  which  would  relieve  m-aconitase  from  inhibition, 
has  been  proposed  as  one  of  the  reasons  for  the  decreased  level 
of  Cit  in  PCa  (30). 

Ch  is  an  essential  component  of  cell  membrane  synthesis  and 
phospholipid  metabolism,  and  functions  as  an  important  methyl 
donor.  Ch-containing  molecules  are  an  essential  component  of 
cell  membranes,  which  are  more  highly  concentrated  in  tumor¬ 
ous  areas  within  the  prostate  than  in  healthy  prostate  tissue 
(32,33).  Ch  groups  are  precursors  and  breakdown  products  of 
the  phospholipid  phosphatidylcholine,  a  major  cell  membrane 


compound  (34).  Increased  Ch  is  observed  as  a  result  of  altered 
phospholipid  metabolism  in  PCa  cell  lines  (33).  This  alteration  is 
most  probably  a  result  of  an  increased  expression  and  activity 
of  choline  kinase,  a  higher  rate  of  Ch  transport  and  an  increased 
phospholipase  activity  (34). 

The  polyamines  Spm,  spermidine  and  putrescine  are  essential 
for  the  differentiation  and  proliferation  of  cells,  the  synthesis  of 
DNA,  RNA  and  proteins,  and  the  stabilization  of  cell  membranes 
and  cytoskeletal  structures  (35).  Previous  studies  have  observed 
high  levels  of  Spm  in  healthy  prostate  tissue  and  benign  pros¬ 
tatic  hyperplasia,  and  reduced  Spm  levels  in  malignant  prostate 
tissue  (16,36-38). 

The  osmoregulator  ml  is  expressed  in  a  variety  of  tissues,  and 
its  decrease  was  observed  in  PCa  within  human  expressed  pros¬ 
tatic  secretions  using  high-resolution  NMR  (39)  and  in  breast 
tumors  (40).  In  our  study,  slightly  increased  ml  ratios  were 
observed  in  cancerous  locations,  but  were  not  statistically 
significant. 

Glu  and  Gin  are  difficult  to  resolve  owing  to  resonance  over¬ 
lap.  As  a  result,  most  MRS  studies  use  the  sum  of  Glu  and  Gin 
(expressed  as  Glx  or  Glu  +  Gin).  Glu  is  extensively  involved  in 
metabolic  and  oncogenic  pathways.  Koochekpour  (41)  showed 
that  serum  Glu  levels  correlated  directly  with  Gleason  scores 
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Figure  5.  Correlation  of  maximum  entropy  (MaxEnt)  and  total  variation  (TV)  non-linear  reconstruction  methods  for  citrate  (Cit)  (A),  glutamine  +  glu¬ 
tamate  (Glx)  (B),  myo-inositol  (ml)  (C)  and  spermine  (Spm)  (D)  in  cancerous  locations,  ft2  values  are  shown  for  each  metabolite. 


Table  1.  Measures  of  sensitivity,  specificity,  positive  predictive  value  (PPV),  negative  predictive  value  (NPV)  and  accuracy  of 
maximum  entropy  (MaxEnt)  and  total  variation  (TV)  methods  using  receiver  operating  characteristic  (ROC)  curve  analysis 


MaxEnt 


Metabolite 

Sensitivity  (%) 

Specificity  (%) 

Accuracy  (%) 

PPV  (%) 

NPV  (%) 

AUC  (%) 

Cit 

86.4 

90.0 

88.6 

90.5 

87.0 

94.0 

Glx 

66.6 

31.8 

47.7 

48.3 

46.7 

56.9 

ml 

54.5 

63.6 

59.1 

60.0 

58.3 

65.1 

Spm 

50.0 

72.7 

61.4 

TV 

64.7 

59.3 

64.9 

Cit 

86.4 

90.9 

88.6 

90.5 

87.0 

92.1 

Glx 

63.6 

36.4 

50.0 

50.0 

50.0 

57.2 

ml 

40.9 

59.1 

50.0 

50.0 

50.0 

58.5 

Spm 

50.0 

63.6 

56.8 

57.9 

56.0 

58.5 

Cit,  citrate;  Glx,  glutamine  +  glutamate;  ml,  myo-inositol;  Spm,  spermine. 


(<6  versus  >8)  and  primary  PCa  aggressiveness.  In  our  study, 
decreased  Glx  ratios  were  found  in  cancerous  locations,  but  were 
not  statistically  significant. 

There  was  an  overlap  between  cancerous  and  non-cancerous 
locations,  possibly  as  a  result  of  the  low  SNR  of  the  Cr  peak.  As 


a  result  of  patient  movement  and  B0  inhomogeneity,  the  reso¬ 
nances  of  Ch,  Cr  and  Spm  are  difficult  to  resolve,  especially  in 
cancerous  locations,  adding  to  the  uncertainty  in  quantification. 
Flowever,  the  use  of  prior  knowledge  fitting  (ProFit)  may  improve 
accurate  metabolite  (42)  quantification,  which  warrants  future 
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investigation.  As  a  result  of  the  limited  patient  population,  we 
did  not  find  any  significant  changes  in  other  metabolites.  In 
addition,  this  study  focused  on  the  PZ  of  the  prostate,  where 
only  80%  of  cancer  occurs.  The  advantage  of  the  compressed 
sensing-based  4D  EP-JRESI  sequence  is  that  it  records  short 
TE-based  spectra  from  multiple  regions  of  human  prostate,  and 
additional  metabolites,  such  as  ml,  Spm  and  Glx,  to  the  normally 
detected  Cit,  Cr  and  Ch. 

This  pilot  work  demonstrated  the  use  of  slice-based  4D 
EP-JRESI,  and  future  work  will  focus  on  volume-based,  five¬ 
dimensional  (5D)  EP-JRESI  in  PCa  (43).  Compressed  sensing- 
based  4D  EP-JRESI  shortens  the  total  acquisition  time,  effectively 
enabling  future  potential  to  extend  to  pathological  evaluations 
in  a  clinical  set-up.  The  current  validation  method  may  require 
further  optimization  to  improve  the  overall  performance.  As 
reported  by  Burns  ef  al.  (26),  the  sample  mask  is  crucial  to  the 
SNR  of  each  reconstructed  prostate  metabolite  for  the  4D  EP- 
JRESI  data.  Further  optimization  of  the  reduction  in  non-linearity 
of  the  reconstructed  peaks  may  enable  accurate  quantification 
of  metabolites.  In  addition,  future  work  will  address  the  use  of 
Poisson  gap  versus  deterministic  sample  masks,  and  the  optimi¬ 
zation  of  the  modulation  functions  for  specific  metabolites 
relevant  to  PCa. 

CONCLUSION 

We  were  able  to  detect  metabolites  in  PCa  using  compressed 
sensing-based  4D  EP-JRESI  data  acquired  in  clinically  acceptable 
times  (<12  min).  We  have  shown  that  it  is  possible  to 
undersample  the  4D  EP-JRESI  sequence  with  an  acceleration  fac¬ 
tor  of  four  times,  and  that  the  data  can  be  reliably  reconstructed 
using  the  TV  and  MaxEnt  methods.  Both  non-linear  reconstruc¬ 
tion  methods  provided  comparable  results. 
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Zohaib  Iqbal1,  Neil  E.  Wilson1,  Rajakumar  Nagarajan1,  Daniel  A.  Margolis1,  Robert  E.  Reiter2,  Steven  S.  Raman1,  and  Michael  Albert  Thomas1 
1  Radiological  Sciences,  University  of  California  -  Los  Angeles,  Los  Angeles,  California,  United  States,  2  Urology,  University  of  California  -  Los  Angeles,  Los  Angeles, 

California,  United  States 

Target  Audience:  Basic  scientists  interested  in  Compressed  Sensing  (CS)  reconstruction  of  semi-laser  based  multi-dimensional  spectroscopic  imaging  (3  spatial  and  2 
spectral  dimensions),  as  well  as  clinicians  interested  in  spectroscopic  imaging  of  in  vivo  prostate. 

Purpose:  Alongside  lung  cancer,  prostate  cancer  is  one  of  the  leading  causes  of  cancer  death  in  men  [1].  J-resolved  spectroscopy  (JPRESS)  is  a  technique  that  utilizes  a 
two-dimensional  (2D)  acquisition  to  spread  overlapping  resonances  into  a  second  spectral  dimension,  and  has  shown  to  be  useful  in  assessing  severity  of  prostate  cancer 
[2,  3].  Spectroscopic  imaging  of  this  technique  provides  greater  spatial  coverage  [4],  however  in  order  to  acquire  data  in  a  clinically  feasible  time,  an  echo-planar 
readout  becomes  a  necessity  to  acquire  one  spatial  and  one  spectral  domain  simultaneously  [5].  In  order  to  further  reduce  scan  time,  a  non-uniform  undersampling 
(NUS)  scheme  can  be  applied  to  the  incremented  dimensions,  and  the  resulting  data  can  be  reconstructed  using  CS  reconstruction  [6].  Due  to  the  limited  bandwidth  of 
classical  180°  pulses,  chemical  shift  displacement  error  (CSDE)  is  apparent.  Using  a  pair  of  adiabatic  180°  pulses  to  replace  these  smaller  bandwidth  pulses  has  shown 
to  limit  CSDE  [7].  The  purpose  of  this  study  was  to  develop  a  5  dimensional  (5D)  technique  (3  spatial  and  2  spectral  dimensions)  capable  of  covering  the  entire  prostate 
utilizing  an  echo-planar  readout,  a  NUS  scheme  on  the  (ky,kz,ti)  dimensions  with  CS  reconstruction,  and  semi-Laser  pulses  to  reduce  CSDE.  This  sequence,  called  semi- 
Laser  5D  echo-planar  J-resolved  Spectroscopic  Imaging  (5D  EP-JRESI),  was  used  to  investigate  the  prostate  for  healthy  volunteers  and  prostate  cancer  patients. 
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Methods:  A  semi-Laser  5D  EP-JRESI  sequence  with  a  maximum  echo 
sampling  scheme  was  developed  and  used  to  scan  four  healthy 
volunteers  (mean  age  =  54.5  years  old)  and  two  prostate  cancer  patients 
(mean  age=  54  years  old)  using  a  multi-channel  body  array  coil  for 
healthy  or  a  single  channel  endorectal  coil  for  patients.  For  the  semi¬ 
laser  based  sequence,  the  overall  structure  was  similar  to  the  5D  EP- 
JRESI  sequence  [8],  however  each  180°  pulse  was  replaced  with  a  pair 
of  180°  adiabatic  pulses  [7].  The  ti  increment  was  placed  in  between  the 
two  pairs  of  adiabatic  pulses.  The  pulse  sequence  diagram  can  be  seen 
in  Figure  1.  The  following  parameters  were  used  for  data  acquisition: 

TE=41  ms,  TR=1000  ms,  spectral  band width=  1190/1000  Hz,  number  of 
t2  points  =  256,  number  of  ti  points  =  64,  FOV  =  16x16x12  cm3,  and 
resolution  =  lxlxl. 5  cm3,  and  the  total  scan  time  was  approximately  17 
minutes.  The  echo  planar  readout  was  used  to  simultaneously  acquire 
32  (oversampled)  kx  points  and  256  t2  points,  and  the  other  spatial 
dimensions  were  acquired  using  16  points  for  ky  and  8  points  for  kz.  A 
NUS  scheme  was  applied  to  achieve  an  acceleration  factor  of  8x.  Before  the  data  were  reconstructed,  the  data  were  truncated  to  only  include  the  spectral  region  of 
interest  (F2  =  1. 2-4.3  ppm)  in  order  to  ensure  minimal  water  contamination.  The  data  was  reconstructed  using  a  modified  split  Bregman  algorithm  that  minimizes  total 
variation  (TV)  [6].  Select  voxels  from  healthy  and  unhealthy  regions  in  the  prostate,  as  determined  from  biopsy,  were  qualitatively  assessed. 
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Figure  1.  The  pulse  sequence  diagram  for  semi-Laser  5D  NUS  Echo-Planar  J-resolved 
Spectroscopic  Imaging  (5D  EP-JRESI). 
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Figure  2.  MRI  localization  (left)  of  axial  (top)  and  coronal  (bottom)  views  of  a 
cancer  patient’s  (age  =  51  years  old)  prostate  are  shown.  A  citrate  metabolite  map 
(middle)  shows  signal  from  three  different  slices.  Two  spectra  (right)  are  shown 
corresponding  to  the  colors  on  the  citrate  metabolite  map.  Red  is  the  healthy  region, 
determined  from  the  biopsy,  and  yellow  is  the  unhealthy  region. 
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Results:  Qualitative  analysis  of  the  2D  spectra  from  the  5 1  year  old  cancer 
patient,  displayed  on  the  right  in  Figure  2,  shows  that  lower  citrate  and 
higher  choline  are  present  in  the  voxel  that  corresponds  to  the  unhealthy 
region  of  the  prostate  (as  determined  from  a  biopsy). 

Discussion:  The  results  show  that  for  this  patient,  healthy  and  unhealthy 
prostate  regions  could  be  differentiated  based  on  the  2D  spectra.  However, 
due  to  motion  during  the  scan,  several  voxels  contained  severe  lipid 
contamination  (hyper-intense  regions  on  the  citrate  metabolite  map).  This 
contamination  may  be  minimized  by  using  an  inversion  pulse  for  lipid 
suppression  [9],  and  will  be  implemented  in  the  future. 

Conclusion:  Semi-Laser  5D  EP-JRESI  has  been  implemented  in  healthy 
volunteers  and  prostate  cancer  patients,  and  is  capable  of  full  prostate 
coverage.  Future  work  will  focus  on  recruiting  more  healthy  volunteers 
and  patients,  quantifying  the  2D  spectra  using  prior  knowledge  fitting,  and 
further  optimizing  the  sequence  to  avoid  lipid  contamination. 
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TEACHING  POINTS 

Multiparametric  magnetic  resonance  imaging  techniques  are  increasingly  being  used  to  address  bottlenecks  in  prostate  cancer  patient  management.  These  techniques  yield  qualitative,  semi-quantitative  and  fully 
quantitative  biomarkers  that  reflect  on  the  underlying  biological  status  of  a  tumor.  1 .  Teach  the  concept  of  prostate  MRI  and  MRSI.  2.  Demonstrate  the  importance  of  functional  and  anatomic  information  in 
prostate.  3.  Show  clinical  applications  that  benefit  from  multiparametric  MRI.  4.  This  functional  and  anatomical  MRI  and  MRSI  guide  urologists  and  radiotherapists  so  that  they  may  be  used  for  individual 
patient  decision  making. 
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